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1.0  INTRODUCTION 


A  detailed  knowledge  of  the  flow  around  an  airfoil  leading  edge  is 
critical  to  the  understanding  of  the  stall  phenomena.  In  an  early  study 
by  McCullough  and  Gault  (1),  various  types  of  stall  were  discussed  and 
extensive  data  were  presented  for  several  NACA  airfoil  sections.  More 
recently,  Roberts  (2)  described  the  phenomenon  and  proposed  a  semiempirical 
theory  for  the  separation  bubble.  Some  investigators  have  tried  to  cal¬ 
culate  the  separated  region  directly  by  using  inverse  boundary  layer  tech¬ 
nique.  These  attempts  have  been  quite  successful  in  obtaining  correlation 
with  experiment  (see  Refs.  3  and  4).  However,  these  methods  are  valid  only 
for  small  bubbles  (less  than  2%  chord)  and  the  application  to  larger  bubbles 
has  not  been  particularly  successful. 

The  current  study  was  undertaken  to  develop  an  analysis  method  for 
laminar  separation  bubbles  (long  or  short)  on  two-dimensional  airfoil  sec¬ 
tions  at  incidence.  A  viscous/potential  flow  iterative  procedure  was  chosen 
due  to  its  simple  and  efficient  nature.  This  view  has  recently  been  con¬ 
firmed  by  Johnson  et  al .  (5);  i.e.,  that  a  good  inverse  boundary  layer  cal¬ 
culation  method  is  as  accurate  as  the  more  time-consuming  Navier-Stokes 
methods  in  predicting  the  separated  flow. 

In  the  present  work,  Cebeci's  boundary  layer  calculation  method  (6), 
which  is  capable  of  predicting  separated  flow  by  an  inverse  boundary  layer 
calculation  procedure,  is  coupled  with  the  potential  flow  calculation  method 
developed  earlier,  VS2D  (7).  The  boundary  layer  procedure  is  a  finite-dif¬ 
ference  method,  sometimes  referred  to  as  the  "Box  Scheme",  and  uses  the  Cebeci- 
Smith  two-layer,  eddy-viscosity  model  for  turbulence  closure.  The  potential 
flow  is  calculated  by  a  low-order  panel  method  where  each  panel  is  represented 
by  a  constant  potential  surface.  The  laminar  separation  bubble  is  modelled 
in  the  potential  flow  calculation  in  such  a  way  that  it  gives  constant  pres¬ 
sure  along  the  surface  inside  the  bubble. 

The  coupled  calculation  procedure  has  been  applied  to  the  NACA  64A006 
airfoil  and  satisfactory  results  have  been  obtained.  The  details  of  the  solu¬ 
tion  procedure  and  comparisons  with  experiment  are  fully  described  in  this 
report. 
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2.0  POTENTIAL  FLOW  CALCULATION  METHOD 

The  velocity  potential  of  a  two-dimensional  body  in  a  uniform  flow 
(unit  velocity  In  x-direction)  In  the  rectangular  Cartesian  coordinate 
system  is  given  by: 


$(x,y)  a  x  +  <fr(x,y) 


(D 


where  $  is  the  disturbance  velocity  potential.  The  solid  surface  boundary 
condition  gives  3 t/3n  =0;  i.e.. 


3£  _  _  3x 
3n  ~  "  3n 


(2) 


From  Green's  identity,  the  disturbance  velocity  potential  on  the  boundary 
can  be  expressed  as 


s 


y  [*t  SnT  (ln  r^>  '  (tn 


(3) 


where  s  and  t  denote  points  on  the  boundary  and  r$t  represents  the  distance 
between  the  two  points  as  shown  in  Figure  1. 


Figure  1.  Integration  Domain. 
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The  integration  domain  of  Eqn.  (3)  includes  the  thin  sheet  which  extends 
from  the  trailing  edge  to  infinity  in  order  to  account  for  the  discontinuity 
in  potential  at  the  trailing  edge. 

Substitute  Eqn.  (2)  into  (3)  to  obtain 


+s  *  l  fh  4;"1"  ?^>dt  ♦  J  /  ^ 


p  =  1  f  p.  -L- )dt  +  -  f  ~  tn(— )dt 

?s  "  J  *t  3nt'  r$t'  tt  J  3nt 


Cl+C2+C3 


Cl+C2+C3 


Since  the  Kutta  condition  requires  upper  and  lower  surface  velocities 
to  be  equal  at  the  trailing  edge,  the  velocity  potential,  $,  along  the  paths, 
C1  and  Cj,  may  be  written  as 

*U  =  *T.E.U  +  It  *  (t  "  ^.E,* 

T‘E*  (5) 

\  =  ^T.E.L  *  It  ’  ^  ‘  V.E.* 

T.E. 

Thus,  the  contributions  from  the  second  term,  •  (t  -  ty  ^  ),  and  from 

T.E. 

the  second  integral  along  the  paths  and  C-j  are  cancelled  out.  This  leaves 
us  with 

[ ♦»  dt 

*«  *T.E.L  y  3n^  r^)dt  *  1  *T.Z.uf  rj^)dt 
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(6) 


This  is  an  integral  equation  of  the  second  kind  and  the  solution  can  be 
found  easily  using  a  simple  integration  technique. 

In  the  present  method,  the  contour  is  approximated  by  a  number  of 
linear  surface  panels  and  t  assunes  a  constant  value  on  each  element.  With 
this  approximation,  the  contribution  from  each  surface  element  can  be  cal¬ 
culated  analytically  and,  hence,  converts  the  integral  equation,  Eqn.  (6), 
into  a  system  of  linear  equations. 


a . .  * .  =  b. 
*3  1 


(7) 


Having  obtained  <p,  the  velocity  on  the  surface  can  be  readily  obtained  by 
a  simple  differentiation. 


u(t) 


3t  1* 

at  "  at  at 


(8) 


It  is  useful 
be  expressed 
i ncorporated 


to  note  that  the  boundary  layer  displacement  effect,  vhich  may 
as  the  velocity  component  normal  to  the  surface,  can  be  easily 
in  the  calculation  by  modifying  the  boundary  condition,  Eqn.  (2). 


at 

an 


V 


or 


3  v  .  1* 
an  3n 


(9) 


where  v  is  the  transpiration  velocity. 
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2.1  Modelling  of  the  Separation  Bubble 


The  separation  bubble  is  represented  by  a  constant  pressure  region  in 
the  potential  flow  calculation;  i.e., 


C 


P 


00) 


or 


01) 


where  f  is  the  total  potential  and  v  is  the  transpiration  velocity.  However, 
if  this  condition  is  imposed  in  addition  to  Eq.  (7),  a  proper  <p  distribution 
may  not  exist  for  a  given  distribution  of  v,  which  satisfies  both  Eq.  (7)  and 
(11)  simultaneously.  In  other  words,  the  system  is  overdetermined.  The 
correct  ^  which  satisfies  both  equations  can  be  obtained  only  if  v  is  com¬ 
patible  with  4>.  In  the  present  work,  this  is  achieved  in  the  following 
manner. 

Consider  Eq.  (3). 

*,  =  i  /[*t3Tt  (tn  ^]dt  (3) 

c 


where 


9<j) 


00 


First,  solve  for  <J>  for  a  given  distribution  of  v.  Here,  Eq.  (3)  or  (7)  and 
(11)  are  used  for  the  attached  flow  region  and  the  bubble  region,  respectively. 
Once  4>  is  calculated,  the  new  v  distribution  can  be  obtained  from  Eq.  (3). 

Note  that  Eq.  (3)  is  an  integral  equation  of  the  first  kind  for  v.  If  this 
v  is  the  same  as  that  of  the  previous  iteration,  then  <j>  and  v  are  considered 
to  be  compatible.  By  repeating  this  procedure,  the  <J>  and  corresponding  v 
distributions  can  be  obtained. 
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3.0  BOUNDARY  LAYER  CALCULATION  PROCEDURE 

3. 1  Basic  Equations 

The  coordinate  system  employed  in  the  boundary  layer  analysis  is  a 
body-fitted  surface  coordinate  system  (x,y),  i.e. ,  one  coordinate,  x,  is 
parallel  to  the  surface  and  the  other,  y,  is  normal  to  the  surface.  This 
is  a  natural  choice  for  the  solution  of  boundary  layer  equations  on  a  two- 
dimensional  airfoil  as  it  doesn't  exhibit  any  singular  behavior  around  the 
leading  edge.  In  such  a  coordinate  system,  the  first-order  boundary  layer 
equations  are: 

Continuity:  =  0  (12) 


Momentum:  [(1  *  e*>  $  J  (13) 


where  e+  is  dimensionless  eddy  viscosity,  e/\>. 

Equations  (12)  and  (13)  are  subject  to  the  usual  boundary  conditions, 
namely. 


u 

o 

II 

> 

II 

at 

y 

u 

=  . 

e  1 
( 

at 

y 

v 

=  0  » 

Now  dimensionless  variables  are  introduced: 
u  =  u/U^,  v  =  v/U^  • 

x  =  x/L,  y  =  y/L  •  ^R|\ 


=  0 

=  6 


P  =  P/pUm2 
Rl  =  U.L/v 


(14) 


(15) 


where  L  and  U,,  are  reference  length  and  velocity,  respectively. 
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Substituting  these  new  variables  in  Eqs.  (12)  and  (13)  yields 


3u 

3x 


(16) 


+  _L 

dr  3y 


(17) 


where  b  =  1  +  e+. 

Introducing  a  stream  function,  defined  as 


u 


(18) 


Equations  (16)  and  (17),  after  dropping  the  bars  for  convenience,  may  be 
written  as 


(K  F"l*  *  +  F'  dF*  _  ci.  3F 

}  dx  *  3x  F  3x 


(19) 


Here  primes  denote  differentiation  with  respect  to  y. 

Equation  (19)  can  then  be  converted  to  a  system  of  first-order  dif¬ 
ferential  equations  and  can  be  solved  by  using  a  finite -difference  scheme 
with  appropriate  boundary  conditions  and  initial  condition.  The  details  of 
the  solution  procedure  will  be  discussed  in  a  later  section.  For  laminar 
flow  calculations,  the  term  e+  in  Eq.  (13)  has  no  meaning  and  is  given  a 
zero  value. 


3.2  Closure  Model  for  Turbulent  Flow 


A  simple  eddy  viscosity  model,  "zero-equation"  model,  is  used  as  a 
closure  relationship.  In  this  model,  Reynolds  stress  is  given  by 


where  e  is  the  eddy  viscosity.  Here  the  variation  of  e  across  the  boundary 
layer  is  prescribed  in  two  parts  as  proposed  by  Cebeci  and  Smith  (Ref .  9  ): 


Inner  Layer 


5i  =  !<*[“  exp(-y/A>]  |  2||H 


£t  i  eo 


(21a) 


Outer  Layer 


=  0.0168 


x> 


f 


u)  dy 


where  <  =  0.4,  A  =  26-j/u  ,  and  u 

t  r 


(21b) 


The  ♦•ransition  effects  are  taken  into  account  by  multiplying  e  by  the 
i ntermi ttency  factor,  vtr,  based  on  Emmon's  hypothesis  that  the  transition 

phenomenon  in  a  boundary  layer  is  characterized  by  the  intermittent  appearance 
of  turbulent  spots  which  move  downstream  with  the  fluid: 


Y 


tr 


=  l  - 


I 

1200 


0.66  „ 

R  (-5- 

xtr  xtr 


(22) 


3.3  Transition  from  Laminar  to  Turbulent  Boundary  Layer 

Two  transition  criteria  are  employed  in  the  boundary  layer  analysis. 
Granville’s  procedure  (Ref. 10)  is  used  for  attached  flow,  and  transition 
in  the  separated  region  is  predicted  by  Crimi /Reeves '  criterion  (Ref.  11). 


3.3.1  Granville's  Transition  Criterion 

Granville  (Ref. 10)  has  developed  a  procedure  based  on  the  relationship 
between  the  neutral  stability  point  and  the  transition  point.  The  neutral 
stability  point  is  defined  as  a  point  downstream  of  which  small  disturbances 
are  amplified  within  the  boundary  layer  and  ultimately  lead  to  transition. 
Smith  (Ref.  12)  and  others  here  proposed  the  correlation  between  the  in¬ 
stability  curve  and  the  local  pressure  gradient,  \  =■  9*/v  •  dU/dx,  as  follows: 
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X  -  -0.4709  ♦  0.11066  in  R0  -  0.0058591  in:R0  0  <  R0  <  650  (23) 

X  =  0.69412  -  0.23992  in  R0  +  0.0205  in:R0  650  <  RQ  <  10,000  (24) 


If  for  a  given  Re,  X  as  calculated  by  Eqns.  (23)  and  (24)  is  greater  than 
that  determined  by  the  boundary  layer  development,  the  flow  has  passed  from 
a  stable  to  an  unstable  region.  Once  the  flow  passes  into  the  unstable 
region,  the  transition  process  begins,  and  Granville  has  been  able  to  show 
that  a  correlation  similar  to  the  instability  process  can  be  used  to  de¬ 
termine  the  transition  point. 

Granville  formed  an  average  pressure  gradient  parameter,  X,  defined  as 


which  correlated  reasonably  well  with  the  momentum  thickness  Reynolds 
number  at  transition  R  .  This  correlation  is  presented  in  analytical 

Htrans 

form  as  follows: 


Transition  Curves 

X  =  -0.0925  +  7.0  x  10'5  R0  (26) 


for  0  <  R„  <  750; 
9tr  “ 
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and  \  =  i. 59381  -  0.45543  in  IT  +  0.032534  in2  Ra 

for  1100  <  Ra  £  3000. 


(28) 


When  the  A  calculated  by  one  of  the  above  expressions  for  a  given  Rg  is 
greater  than  the  value  determined  from  the  boundary  layer  development, 
transition  is  predicted. 


3.3.2  Crimi /Reeves 1  Criterion 

If  the  laminar  boundary  layer  separates  prior  to  transition,  then  the 
following  relationship  is  used  to  predict  the  onset  of  transition  as  proposed 
by  Crimi  and  Reeves  (Ref.  11), 


y 

u=0  _ 

sep 


(29) 


where  5*  is  the  boundary  layer  displacement  thickness  and  the  subscript,  "sep 
denotes  the  quantity  at  the  separation  point  (see  Figure  2). 


3.4  Finite-Difference  Scheme 


(33) 


i  l 

Ui  -  UJ-l 
hj 


(V1)1  +  (V1)1'1  _  u1  -  u1"1 
2  ~  k< 


(34a) 


Equation  (34a)  can  be  further  discretized  as 


,  /  i  i  i-1  i - 1  \  ,  /  i  l  i-l  i-l\ 

%(vj '  Yi  *  vj  -  vj-i)  "  2k7(uj  *  uj-i  *  uj  -Vi) 


i-l  i-l 


(34b) 


or 


i  i  h.  ,  i  i  v  i-l  i-l  Zh. 

V  Vl  ‘  ( V  Vl  )  ’ ‘'j 


i-l  i-l  2h . 


(34c) 


Equations  (33)  and  (34c)  together  with  boundary  conditions  form  a  block  tri- 
diagonal  matrix  equation  and  can  be  readily  solved  for  U  and  V  by  the  block 
elimination  method  (Ref.  13). 

The  accuracy  of  this  method  is  of  second  order  and  its  implicit  nature 
provides  unconditional  stability  for  any  choice  of  x-step  size. 


3.5  Solution  Procedure 


3.5.1  Standard  Boundary  layer  Calculation 

As  shown  earlier,  the  present  "Box"  scheme  requires  the  differential 
equations  to  be  of  first  order.  By  introducing  new  independent  variables, 
u(x,y)  and  v(x,y),  Eq.  (19)  can  be  transformed  into  a  system  of  first- 
order  equations  as  follows: 
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F*  »  u 


t 


h 


U*  *  V 


(»,»)• 


ft* 


9u 

3x 


(35) 


The  usual  boundary  conditions,  non-slip  condition  at  the  wall  and  free- 
stream  condition  at  the  outer  edge  of  the  boundary  layer,  are  given  by 


F  *  0,  u  =  0 

at 

o 

ii 

>. 

u  =  Ue 

at 

y  =  6 

(36) 

These  represent  basic  equations  for  the  standard  boundary  layer  calculation 
procedure  and  can  be  readily  solved  by  the  method  described  earlier. 


3.5.2  Inverse  Boundary  Layer  Calculation 

Equations  (35)  become  singular  at  the  separation  point  and,  therefore, 
the  above  solution  procedure  fails  as  the  flow  approaches  separation.  This 
difficulty  can  be  avoided  if  a  different  set  of  boundary  calculations  are 
specified;  i.e.. 


F  =  0,  u  =  0  at  y  =*  0 

F  ■  U@(y  -  6*)  at  y  ■+■  ■  (37) 

In  this  inverse  mode,  Ue  comes  out  as  part  of  solutions  and  need  not 
be  specified.  This  is  an  iterative  procedure  that  is  repeated  until  the 
calculated  6*  matches  the  specified  one. 

For  flows  with  negative  wall  shear,  it  is  necessary  to  make  approximations 
to  the  governing  equations  to  continue  the  calculations  past  the  separation 
point.  The  approximation  used  here  neglects  the  u(au/3x)  term  in  the  region 
of  negative  u-velocity  as  originally  suggested  by  Reyhner  and  Flugge-lotz 
(Ref.  14). 


V'- 
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The  calculation  starts  from  the  forward  stagnation  point  with  the 
Hiemenz-Howarth  two-dimensional  stagnation  profile  (Ref.  15)  as  an  initial 
condi tion. 


3.5.4  Grid  Distribution  across  the  Boundary  Layer 


The  computation  domain  in  the  direction  normal  to  the  surface  extends 
from  the  wall,  where  the  no-slip  condition  applies,  to  some  point  beyond  the 
edge  of  the  boundary  layer,  where  the  velocity  component,  u,  approaches  its 
potential  flow  value,  Ue.  Typically,  51  to  61  grid  points  are  distributed 
across  this  domain.  Additional  grid  points  are  added  as  the  boundary  layer 
grows  to  ensure  that  the  whole  boundary  layer  lies  within  the  integration 
domain.  The  node  points  are  redistributed  whenever  the  number  of  points 
reaches  the  prescribed  maximum  (61  in  the  present  case). 


The  distribution  of  these  grid  points  is  nonuniform,  with  the  spacing 
denser  near  the  surface.  The  use  of  nonuniform  grid  spacing  across  the  layer 
is  essential  for  turbulent  flow  calculations.  This  can  be  achieved  by 
using  a  grid  which  has  the  property  that  the  ratio  of  lengths  of  any  two 
adjacent  intervals  is  a  constant;  i.e.,  h.  *  Ch^.  The  distance  to  the 

jth  line  is  given  by  the  following  formula: 


=  hj(C  -  1 ) / (C  -  1) 


C  >  1 


The  total  number  of  points  J  can  be  calculated  from  the  expression 


ZnC 


where  yj  is  the  outer  limit  of  the  integration  domain. 
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4.0  VISCOUS/POTENTIAL  FLOW  ITERATION 


Initially,  the  potential  flow  solution  is  obtained  with  or  without 
the  separation  bubble  specified  by  the  procedure  described  in  an  earlier 
section.  This  potential  flow  solution,  mainly  velocity  distribution,  is 
used  in  the  calculation  of  the  boundary  layer  development.  The  attached 
flow  solution  starts  from  the  stagnation  point  and  proceeds  until  it  predicts 
the  separation.  If  a  turbulent  separation  is  predicted,  then  the  boundary 
layer  calculation  stops  and  the  calculation  returns  to  the  potential  flow 
routines.  If  the  calculation  predicts  a  laminar  separation  along  the  upper 
surface  of  the  wing,  the  solution  continues  with  the  inverse  boundary  layer 
method.  The  boundary  layer  displacement  thickness  distribution  which  is 
required  for  the  inverse  calculation  is  prescribed  i ni ti ally— later  this  will 
be  correlated  with  the  normal  velocity  distribution  in  the  potential  flow 
calculation— so  that  the  calculation  can  continue  until  the  flow  reattaches. 
Reattachment  is  possible  because  of  transition  to  the  turbulent  flow  inside 
the  bubble.  The  section  downstream  of  the  reattachment  point  is  calculated 
by  a  standard  boundary  layer  method.  This  completes  one  full  viscid/i nviscid 
iteration  cycle  (see  Figure  4). 

The  calculation  returns  to  the  beginning  and  repeats  the  procedure.  Now 
the  normal  velocity  due  to  the  boundary  layer  displacement  effect  can  be  in¬ 
corporated  in  the  potential  flow  calculation.  The  whole  procedure  is  con¬ 
sidered  to  be  converged  when  the  external  velocity  along  the  bubble  matches 
that  of  the  potential  flow  and  the  reattachment  point  does  not  move  from  one 
iteration  to  another. 


Figure  4.  Schernetic  of  Cslculfltion  Procedure. 
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5.0  RESULTS  AND  DISCUSSION 


The  experimental  study  of  laminar  separation  bubbles  on  an  airfoil 
section  has  been  pursued  by  only  a  small  number  of  investigators  and,  con¬ 
sequently,  few  adequate  data  sets  for  comparison  are  available  in  the 
literature. 

Gaster  (16)  generated  a  separation  bubble  on  a  flat  plate  by  artificially 
applying  an  adverse  pressure  gradient  and  was  successful  in  measuring  the 
characteristics  of  laminar  separation  bubbles;  i.e.,  pressure,  velocity  pro¬ 
files,  etc.  However,  since  one  of  the  prime  objectives  of  the  present  study 
is  modeling  of  the  separation  bubble  in  the  potential  flow  calculation,  this 
case  was  not  considered.  The  data  by  Gault  (17)  for  various  airfoil  sections 
are  limited  to  small  bubbles  and  not  suitable  for  the  present  work.  In  1951, 
McCullough  and  Gault  conducted  an  experimental  study  on  the  leading-edge  bubble 
for  various  airfoil  sections.  Among  the  data  obtained  by  McCullough/Gault, 
the  data  on  the  airfoil  section  NACA  64A006  were  chosen  for  comparison  because 
of  the  medium  range  bubble  sizes  (3  -  20%)  which  are  of  current  interest. 

Figures  5  and  6  are  the  results  of  two  full  viscous/potential  flow 
iterations  for  a  =  5°  and  6°,  respectively.  Figure  7  shows  the  pressure  dis¬ 
tribution  for  a  =  7°  just  after  the  potential  flow  calculation  with  the  sepa¬ 
ration  model.  All  the  solutions  are  started  with  the  flow  fully  attached  and 
the  boundary  layer  analysis  (both  direct  and  inverse  methods)  determines  the 
laminar  separation,  transition  and  reattachment  points.  In  this  separation 
bubble  analysis,  the  location  of  transition  is  very  important  as  the  bubble 
model  is  applied  from  the  separation  point  to  the  transition  point.  The 
transition  procedure  adopted  here  (Crimi /Reeves)  tends  to  predict  early  transi¬ 
tion.  It  is  possible  that  the  velocity  profile  on  which  the  separation  cri¬ 
terion  is  based  may  be  in  error  or  the  procedure  itself  may  not  be  adequate 
for  long  bubbles.  This  subject  could  not  be  examined  thoroughly  during  the 
current  phase  of  the  work  due  to  limited  resources  and  deserves  further  study. 

In  order  to  avoid  any  uncertainties  in  this  repect,  the  onset  of  transition 
is  fixed  in  the  present  calculations. 

As  shown  in  Figures  5  and  6  the  boundary  layer  solution  for  the  first 
iteration  is  not  close  to  either  the  potential  flow  solution  on  the  experimen¬ 
tal  data.  However,  with  the  aid  of  the  separation  model,  these  solutions 
quickly  converge  in  the  second  iteration.  Figure  7  shows  a  slight  discrepancy 
in  the  transition/reattachment  region  between  the  potential  flow  solution  and  the 
data.  This  is  due  to  the  transition  process  which  is  not  an  abrupt  process  as 
modelled  in  the  potential  flow  calculation.  This  view  is  confirmed  by  the 
fact  that  the  boundary  layer  solutions  are  in  better  agreement  with  the  ex- 
perimenta  data,  as  a  result  of  the  more  gradual  transition  process,  than  is 
the  potential  flow  calculation  with  its  separation  model. 

As  pointed  out  earlier,  it  is  apparent  from  these  results  that  the  size 
of  the  bubble  is  of  critical  importance  to  the  prediction  of  the  correct 
pressure  level  along  the  surface  of  the  separation  bubble.  It  is  also  use¬ 
ful  to  note  that  the  flow  does  not  exhibit  trail ing-edge  separation,  and  the 
pressure  correlation  for  the  rest  of  the  airfoil  is  very  good  as  shown  in 
Figures  8  through  10.  The  C„  -  a  curve  comparison  is  plotted  in  Figure  11  and 
the  correlation  is  very  encouranging. 
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Figure  6.  Pressure  Distribution  after  Two  Full  Viscous/Potential 
Iterations;  NACA  64A006,  a  =  6°,  Re  =  5.8  x  10°  (Fixed 
Trans i tion) . 
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Figure  7.  Effect  of  Separation  Model  in  the  Potential  Flow 
Calculation. 
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Figure  10.  Overall  Pressure  Distribution  on  NACA  64A006; 
a  *  7°. 
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11.  Ci  -  a  Curve  for  an  Airfoil  Secti 
NACA  64A006 ;  Re  =  5.8  x  106. 


The  overall  results  presented  here  are  quite  good.  Even  better  corre¬ 
lation  might  have  been  obtained  had  the  calculation  gone  one  more  iteration. 

One  critical  area  that  needs  further  study  is  the  criterion  for  identifying 
the  onset  of  transition  to  turbulent  flow  inside  the  bubble.  Relying  completely 
on  the  displacement  thickness  at  the  separation  point  and  the  local  velocity 
profile  as  is  done  by  Crimi /Reeves  may  make  the  overall  calculation  unstable 
as  the  resulting  velocity  profiles  in  the  inverse  boundary  layer  calculation  are 
very  sensitive  to  the  boundary  layer  displacement  thickness  distribution. 

Perhaps  the  boundary  layer  displacement  prescribed  along  the  bubble  should  be 
more  closely  related  to  the  convergence  criteria.  The  effect  of  the  boundary 
layer  displacement  thickness  distribution  on  the  velocity  profile  in  the  separa¬ 
tion  zone  needs  to  be  examined  more  carefully  in  a  future  study. 


6.0  RECOMMENDATIONS  FOR  FUTURE  WORK 


It  has  been  demonstrated  in  this  work  that  the  leading-edge  separation 
bubble  can  be  analysed  by  a  viscous/potential  flow  iterative  procedure.  As 
described  in  the  previous  section  the  correlation  with  the  data  is  encouraging. 
However,  there  is  considerable  room  for  improvement  in  the  calculation  procedure 
to  make  it  more  useful  and  practical. 

The  future  success  of  the  calculation  procedure  hinges  on  the  accurate 
prediction  of  the  onset  of  transition.  The  transition  criterion  inside  the 
bubble  needs  to  be  improved  as  the  current  one,  which  is  based  on  the  relation¬ 
ship  between  the  local  velocity  profile  and  the  boundary  layer  displacement 
thickness  at  the  separation  point,  predicts  early  transition.  There  are  two 
possible  reasons  for  this  premature  transition.  Firstly,  the  Crv  'Reeves  pro¬ 
cedure  itself  may  not  be  valid  for  the  long  bubble.  Secondly,  the  velocity 
profile  obtained  from  the  inverse  boundary  layer  calculation  in  the  region  of 
interest  may  be  in  appreciable  error  and  thus  invalidates  the  transition  test. 
These  uncertainties  can  be  clarified  by  the  following  means:  a  correlation  study 
between  the  Crimi /Reeves  criterion  and  the  existing  long  bubble  data,  and  a 
sensitivity  test  on  the  boundary  layer  velocity  profiles  against  the  displace¬ 
ment  thickness  distribution.  It  may  also  be  desirable  to  compare  the  present 
boundary  layer  method  with  other  methods,  e.g.,  Horton's  inverse  boundary  layer 
method  (18). 

The  speed  of  convergence  of  the  current  iterative  procedure  in  obtaining 
the  potential  flow  calculation  is  relatively  slow.  More  work  is  needed  in 
this  area  to  optimize  the  convergence  process. 
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